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Abstract 

We use a variational method based on the principle of least ac- 
tion to obtain approximate time-dependent single soliton solutions to 
the KdV equation. A class of trial variational functions of the form 
u{x, t) = —A{t) exp —(3{t) \x — q{t)\ " , with n a continuous real vari- 
able, is used to parametrize time- dependent solutions. We find that 
this class of trial functions leads to soliton-like solutions for all n, mov- 
ing with fixed shape and constant velocity, and with energy and mass 
conserved. Minimizing the energy of the soliton with respect to the 



parameter n, we obtain a variational solution that gives an extremely 
accurate approximation to the exact solution. 

1 Introduction 

In a previous work |T]|, we introduced a post-Gaussian variational approxi- 
mation, a continuous family of trial variational functions more general than 
Gaussians, which can still be treated analytically. We demonstrated their 
usefulness by discussing several features of the general nonlinear Schrodinger 
equation: we derived an approximation to the one-dimensional soliton so- 
lution and showed the "universality" of the critical exponent for blowup in 
the supercritical case; for the critical case we calculated the critical mass 
necessary for blowup. 

In the present work we consider another application of the post-Gaussian 
variational method, deriving an approximation to the one-dimensional KdV 
soliton. Of course the soliton solution to this exactly integrable system is 
well-known and easy to find by direct integration. The aim here is rather 
to provide another illustration of the power, generality, and accuracy of our 
method in a problem with some new features, in preparation for treating 
other more difficult problems. 

The time- dependent variational principle is very useful for quantum 
mechanics [Q and field theory problems as well as for studying nonlinear sys- 
tems described by a Lagrangian ||^. In our recent work we showed how one 
could obtain information about solitons and self- focusing (blow-up) |I[], [Q 
in the nonlinear Schrodinger equation using a version of Dirac's variational 
principle and a canonical set of variational parameters. This choice of param- 
eters led automatically to a conserved energy for the effective Lagrangian. 
For the KdV equation, however, a similar choice of variational trial functions 
leads to a Lagrangian that is not automatically in canonical form. We em- 
ploy the method of Faddeev and Jackiw f^ to explicitly obtain a canonical 
Lagrangian with symplectic structure, which leads to a conserved energy for 
the variational parameters as well as a Hamiltonian structure. Assuming 
a variational ansatz of the form u{x,t) = —A{t)exp —Pit) \x — q{t)\ " ,we 
find that the conserved energy of the effective Lagrangian is minimized for 
n = 0.877, giving an extremely accurate energy, velocity, and shape for the 
approximate single soliton solution. 



2 Action for the KdV Equation 
and the Variational Principle 

The starting point for the derivation of the KdV equation is the action 

T = f Ldt, (1) 

where L is given by P 

L = J y-^^x^t - i^xf - ^i^xxfj dx. (2) 

From this Lagrangian, one can determine the conserved Hamiltonian 



H = f [{nif) -L]dx= I 



{'{'xf + -{Vxxf 



dx, (3) 



where vr = ^. Because the action is stationary with respect to variation in 
if (i.e. 5V /5ip = 0), we obtain the equation 

^xt - Q^xVxx + ^xxxx = 0, (4) 

which we recognize as the KdV equation with the identification u = ip^: 

Ut - 6uu^ + u^^x = 0, (5) 

r is the action of the system described by ip{x, t); the variational principle 
is a version of Hamilton's least-action principle. In the variational principle (p 
is an arbitrary square integrable function. To obtain an approximate solution 
to the KdV equation, we consider a restricted class (pv{x, t), constrained to a 
form that tries to capture the known behavior of the full if for the problem 
at hand. In this article we will consider only trial wave functions that are 
capable of describing the motion of an initial configuration which can be 
reasonably approximated at time t = by ^(a;, 0) = — Aexp [— /3|a; — ap"] 
for suitable choice of A, /3, a and n. Thus, we choose for our trial wave 
function 

M„(x, t) = -A{t) exp [-/3{t) \x - qit)M , (6) 

where n is an arbitrary continuous, real parameter. 



The variational parameters have a simple interpretation in terms of ex- 
pectation values with respect to the "probability" 



[u^{x,t)y 



M{t) 



P{x,t) 

where the mass M is defined as 

M(t) = I [u^{x,t)fdx. 



(7) 



(8) 

(Here we allow M to be a function of t, even though M is conserved for the 
exact KdV equation.) 

Since {x — q{t)) = 0, q{t) = (x). From (8) we have 



A{t) 



^l/2(2;5)l/4n 



2rfe + i 



1/2- 



(9) 



The inverse width f3 is related to 



4n/3 



Thus, we can write 



U = ifx 



MV2(2/3)V4n 



exp 



-P{t)\x-q{t)\ 



2n 



(10) 



(11) 



2r(2k + i); 

To evaluate the action we first must determine (p. A convenient choice of 
integration constant gives 



V9(x,t) = -A{t) exp[-/3(t)|y-g(t)|Hrfy 

Jq{t) L J 

= -A(t)(2n)-V'/'"e(a; - q{t))F [/^^/^nj^ _ ^^^^n ^ (^2) 



where 



/5^/2"(x-g(t)) 



I3\x-q{t)[ 



e ^2;2n ^dz. 



(13) 



The sign function e{x) explicitly displays the oddness properties of ip: 



^{x- qit),t) = -(fi (- [x - q{t)] , t) . 



(14) 



For n = 1 we have 



/3^/2 (a; _ g(i))l = ^1/2 Erf [/3V2 (^ _ g(^)) 



(15) 



The oddness property of ip guarantees conservation of M and (3, as we will 
see below. (For a two soliton ansatz neither M nor (3 would be time- inde- 
pendent.) 

We next evaluate the action, given by (|1]) and (|^), for the trial wave 
function in (M). First consider 



'" = -«^"+(lf + 4^)»' + '^''<^-''™' 



where 



(16) 



rx — q(t) 

G{x-q{t)) = A l^r^exp 

Jo 



-m\yr dy. 



(17) 



We notice that G is also odd in x — q{t). Thus, in determining / y^ip^ipA dx, 
only the first term of ([T6|) contributes and we obtain —Mq/2 for the resulting 
integral. This oddness property causes /3 to be a constraint rather than a 
dynamical variable in the ensuing Lagrangian dynamics. 
Evaluating the other terms in (0), we obtain 



r(g,/?,M,n) 



where 



— Mq + Ci(n)/5i/^"M3/2 - C2(n)M/3^/" ) dt 



= jLi{q,q,M,(3)dt, 



(18) 



Ci(n) 



C2(n) 



l/4n 



2r — + 1 

' 2n 



4 r(^ + i) 



-1/2 



(19) 



The variational equations 5T /5qi = yield a set of equations for the 
parameters q, (3 and M. However it is more useful to put our first Lagrangian 
Li into canonical form, so that the conserved Hamiltonian can be displayed 
and the unconstrained Lagrangian can be obtained. Following the work of 
Fadeev and Jackiw |^, we first rewrite Li as 



\(,M 



qM)-H2{(3,M), 



where 



H2{f3,M) = -Ci(n)/5^/^"M3/2 ^ c'2(n)M/3^/". 



(20) 



(21) 



We next recognize that /3 is a variable of constraint since f3 does not appear 
in L2. We eliminate /3 (using 6T/6P = 0) and find 



/3=[(i(n)f"M2"/3, 

where d{n) = [Ci{n)/AC2{n)Y^^. 

We now eliminate (3 in (^) in favor of M to obtain: 



(22) 



where 



L3 = \{qM-qM)-H3{M), 



HsiM) = {C2d^ - Cid) M^/\ 



(23) 



(24) 



We now have unconstrained Lagrangian dynamics with a conserved Hamil- 
tonian H3, which is the Hamiltonian (^ evaluated with the trial wave func- 
tion. The set of Lagrange equations resulting from (|23|) are 



M = =^ M = const. 
=^ (3 = const. 



and 



^ 3 



as well as a conserved energy 



Cid - C2rf^ 



(25) 



(26) 



E = (Carf^ - Cid) M^/^ (27) 

Thus the velocity of the sohton is constant and related to the conserved 
energy via 



Thus we have 



->'- 



(28) 



u,ix,t) = -rf(n)M2/32i/4" 



2r 



i + — 

2n 



-1/2 



X 



expJ-d^'^M^'^/^i^-ct- 



Xq 



\2n 



(29) 



with c given by (^) and (pS]). To find the "best soliton" for our class of trial 

wave functions we need to minimize the energy of the soliton with respect 

to n. In Fig.l we plot the dimensionless energy E{n)/M^^^ . We see that 

the energy is minimized for n = 0.877; at that value we find that E = 
-0.3925M5/3_ 

The exact single soliton solution of the KdV equation is given by H 



u{x,t) 



-sech^ 



f-c^/^(x-ct-Xo)j 



(30) 



2 V2 

where c = {3/2f/^M^/^ = I.SIOM^/^. Using (|), we find that the energy of 
the exact soliton solution is given by 



E 



Kl) 



5/3 



M^/^ = -0.3931M^/^ 



(31) 



Comparing ( |27D and (0), we find that the energy of the true soliton is lower 
than that of the variational approximation (as it must be), but the error is 
only 0.15%. Using the proportionality of c and E, we find that the variational 
calculation gives for the soliton velocity at the optimal n: 

vo = 1.308M2/3 (32) 

which is accurate to 0.2%. 

In Fig. 2 we plot the best variational soliton wave function ul with 
n = 0.877 and the exact single soliton solution u^. To show the accuracy in 



7 



detail we plot ul{x) -~u'^{x) in Fig. 3; the worst error (at the origin) is about 
0.5%. 

From these calculations we conclude that by optimizing the class of non- 
Gaussian variational vave functions we get an extremely accurate approxi- 
mate soliton with correct velocity and energy. 
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